Generic nodeless Larkin Ovchinnikov states due to singlet-triplet mixing 
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Larkin- Ovchinnikov (LO) states typicaUy have a singlet-gap that vanishes along real-space lines. 
These real-space nodes lead to Andreev midgap states which can serve as a signature of LO pairing. 
We show that at these nodes, an odd-parity, spin-triplet component is always induced, leading to a 
nodeless LO phase. We find the two-dimensional weak coupling, clean limit s-wave phase diagram 
when this spin-triplet part is included. The triplet component is large and increases the stability of 
the FFLO phase. We also show that the spin-triplet contribution pushes the midgap states away from 
zero energy. Finally, we show how our results can be explained phenomenologically though Lifshitz 
invariants. These invariants provide a simple approach to understand the role of unconventional 
pairing states, spin-orbit coupling, and inhomogeneous mixed singlet-triplet states that are not due 
to a FFLO instability. We discuss our results in the context of organic superconductors. 



There are strong reasons to suspect that the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) P, 0] 
phases appear in the quasi-one-diniensional(QlD) 
Bechgaard salts (TMTSF)2X[i, i] and in the 
quasi- two-dimensional(Q2D) organics k-(BEDT- 
TTF)2Cu(NCS)2[i and A-(BETS)2GaCl4f65. FFLO 
phases have also been argued to be of importance in 
understanding ultracold atomic Fermi gases 0, Q and in 
the formation of color superconductivity in high density 
quark matter Q. The understanding of these phases 
has become a relevant and topical pursuit in physics. 
A central result of theoretical studies is the ubiquitous 
appearance of the LO phase, a striped superconducting 
phase in which the s pin - singlet order parameter vanishes 
spatially along lines jl(]|. Indeed, it has been suggested 
that the observation of Andreev bound states localized 
at these nodes would provide strong evidence for LO 



phase[ll|,ll2| 



Here we argue that the spin-singlet LO phase is gener- 
ically nodeless due to the appearance of a spin-triplet 
component at the spatial nodes of the spin-singlet com- 
ponent. We further show that the triplet component is 
stabilized by "removing" the Andreev bound states, that 
is, by pushing these states away from zero energy. 

We begin with a microscopic derivation of our main 
results. This derivation considers a 2D superconductor 
with spin-singlet s-wave and spin-triplet p-wave pairing 
interactions. This is followed by a phenomenological de- 
scription that shows how Lifshitz invariants (LI) account 
for the microscopic results and allow for a significant 
generalization to include the effects of unconventional 
pairing states, spin-orbit coupling (SOC), and inhomo- 
geneous singlet-triplet mixed states not due to an FFLO 
instability. While there have been prior studies of the 
role of p-wave interactions on the FFLO phase 0, [l3- 17 
and in a related phase in cold atoms [18[, these studies 
have focussed on the high field region near the normal to 
superconducting phase transition, where the gap is small. 
Here we examine the low-field transition from a usual su- 
perconductor to a LO phase which requires a solution of 
the non-linear Eilenberger equations. 

We use the Eilenberger equations as presented by 



Alexander [ill[Tl,[li. The central equation for the qua- 
siclassical Green's function g{R, k; ten) is 



[ie„f^ - A - w, g] + ivf • Vg = 0. 



(1) 



where k is the direction of the Fermi momentum, e„ = 
7rT(2n -f 1) are the Matsubara frequencies, v f \s the 
Fermi velocity. We denote the three Pauli matrices in 
particle- hole space by Tx,Ty,Tz, and in spin space by 
(ax,o'y,az)= (T. The Green's function must satisfy Eilen- 
berger's normalization condition — — tt^I. The quasi- 
classical Green's function in Nambu space is 

p,^ ( 9 + (/ + /• (T)i(Jy \ , . 

The Zeeman coupling with the magnetic field is given by 

where is the magnetic moment of the electron. The 
order parameter matrix in Nambu space is 



A(Rk)-( ° + ^ • '^^'"y 

The self-consistency relations are 
A{R,k) = No7rTY^ < V(k,k')f{R,k';ien) >^ 



(4) 



(5) 



A(i2, k) = TVo^T ^ < V{k k')fiR, k'- ien) (6) 

n 

where V{k, k') is the pairing interaction, A^o is the den- 
sity of states at the Fermi level, and <>^' denotes the 
average over the Fermi surface. To determine which 
phase is stable, we use the free energy derived from the 
Luttinger-Ward functional by Voronstov and Sauls [20| : 

Af{R)^\j\\TY,Noj^TrA{g^-\g), (7) 
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g\ is an auxiliary propagator obtained from the solution 
to the Eilenberger equation with the physical order pa- 
rameter scaled by the dimensionless coupling parameter 
< A < 1, 



[i^nTz - AA - v,g\] 



V5A = 0. 



(8) 



We include both singlet s-wave interactions and triplet 
p-wave interactions. We assume a 2D cylindrical Fermi 

surface and a paring interaction V{k, k') = Vs + Vtk ■ k' . 
The relative strength of triplet interaction is given by the 
parameter Tp — Tt/Tg where Tg (Tt) are the for the 
singlet (triplet) pairing. Due to spin rotational invari- 
ance, we will get equivalent results for the field chosen 
along any direction. We therefore set the field along z 
direction for convenience. However, we note that the 
magnetic field should be in the plane to ensure that vor- 
tices can be ignored. Similarly, we also assume spatial 
variations along the x direction. The structure of the 
Eilenberger equations then ensure that there will be a 
non-zero spin triplet component of the order parameter 
of the form d = z^/2kx/kfipz{x) = z^/2cos9iIjz{x). More 
specifically the self-consistency relations become 



2^ 



(9) 



i;z{R) ^ NonTVtY] / V2 
„ Jo 



cos(e)^/. 



■le, 



0- 



(10) 

that appears 



This leads to the gap function -0 -I- azkipz {x) 
in the Eilenberger equations for g. 
Phase Diagram 

In the vicinity of the transition from the normal 
state to the superconducting states, we set (ip^ipz) = 
e^"^^ {ip,ipz) and find the instability line Hc2 by solving 
the linear gap equation and optimizing Hc2 with respect 
to q. The order parameter for a particular q, ipq, is a 
linear combination of the singlet and triplet parts, that 
is {ip,ipz) = {oL, j3)ij}q. Due to parity symmetry, this 
solution has the same Hc2 as V'-gj which is given by 
(■0, ipz) = (a, —l3)ijj~q- As a consequence, just below Hc2, 
two solutions can appear: a solution for which only one 
of ipq or ip^q is non-zero (known as the FF phase); or 
a solution for which both are non-zero and \ipq\ = \tp-q\ 
(known as the LO phase). To determine which of these 
phases appear at Hc2 requires an analysis beyond the 
non-linear gap equation. Keeping up to order in the 
free energy, we find that both the FF and the LO phases 
appear. The FF phase takes up only a small portion 
of the phase diagram. Nevertheless, this has an impor- 
tant physical consequence. In particular, if the FFLO 
phase is generated created by a magnetic field applied 
in the plane, then an additional magnetic field applied 
along the z direction will lead to vortices. The degener- 
acy of the FF and LO phases ensures that there exists a 




FIG. 1. Singlet (circles-solid) and triplet (squares) order pa- 
rameters at r = 0.2Ts for 2D FFLO superconductors with 
Tp = 0.5 and = Vf/{2TvTc). 



stable vortex lattice of half- quantum vortices, as opposed 
to the usual Abrikosov lattice of full-quantum vortices 
2l| . This half-quantum lattice will exist in a region near 
where these two phases are degenerate. 

We also compute the phase boundary from the uniform 
superconducting phase to the FFLO phase. In general, 
this requires a numerical solution of the Eilenberger equa- 
tions. We use an efficient and numerically stable method 
described by SchopohlQ in which the Eilenberger equa- 
tions are transformed to Riccati equations. The transi- 
tion from the the uniform superconducting phase to the 
FFLO phases is found by computing the free energy of 
these two phases. Fig. 1 shows the self-consistent order 
parameter at the transition from the uniform supercon- 
ducting state to FFLO state for Tp = 0.5 and T/Ts = 0.2. 
The spin-singlet order parameter is qual itatively similar 
to previous results on the LO phase[ll[. However, the 
spin-triplet order parameter is maximum where the spin 
singlet order parameter vanishes, removing the spatial 
line nodes usually predicted in the LO phase. Further- 
more, we find that if the spin-singlet order parameter 
is chosen real, then the spin-triplet order parameter is 
imaginary. Both the phase and and the positions of the 
maxima of the spin-triplet order parameter are a natu- 
ral consequence of the phenomenological arguments pre- 
sented later. The complete H-T phase diagrams are pre- 
sented in Fig. 2 for Tp = 0.0 and Tp = 0.5. 



Quasiparticle Properties 

Previous studies of the LO phase have found midgap 
Andreev states associated with sign change of the spin- 
singlet order parameter 11 , l3] . Given the removal of the 
gap through the appearance of a spin-triplet order pa- 
rameter, wc compute the single particle density of states 
to see what happens to these midgap states. The local 
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FIG. 2. FFLO phase diagrams for Tp = Tt/Ts = and Tp = 
Tt/Ts = 0.5. At low fields, the uniform superconducting to 
LO phase transition is second order (circles-dot). When Tp = 
Tt/Ts = 0.5, a FF phase appears in a small region of the 
phase diagram (solid lines). 



quasiparticle density of states (LDOS) at point R with 
spin direction e can be calculated from 

Ne{R; e) = - < -Im{g{R, k;e) + e- g{R, fc; e)) >. 

(11) 

where ie„ — ^ e + iO"'". In Fig. 3, we show the LDOS at the 
nodes of the spin-singlet order parameter for spin-up ex- 
citations. These results compare the solutions for Tp = 
and Tp = 0.5. The LDOS for spin-down electrons can be 
found by reflecting LDOS for spin-up electrons through 
zero energy. When Tp = 0, there exist Andreev bound 
states with energies pinned to the middle of the gap. This 
agrees with previous studies 



Il|,|l2|. Once the spin-triplet 




FIG. 3. Spin-dependent local density of states (LDOS) at 
nodes of spin-singlet order parameter for Tp = (dot-dashed) 
and Tp = 0.5 (solid lines, the y axis has been offset by 1 for 
clarity) . 



zero energy. This shifting of these provides a microscopic 
mechanism through which the spin-triplet order param- 
eter is energetically stabilized. We note that a similar 
Andreev bound state removal mechanism has been pro- 
posed to explain the occasional appearance of spin den- 



sity wave (SDW) order at the spin-singlet nodes [2 3l 
An important difference with our results is that the spin- 
triplet order we find is required to appear by symmetry 
while the SDW order is not. 

One physical property associated with the Andreev 
midgap states is the appearance of an increased ferromag- 
netic magnetization at the nodes of the spin-singlet order 
parameter . To investigate the role of the spin-triplet 
order parameter on this, we calculate the magnetization 
and find that the spatial peak of magnetization and the 
total magnetization both decrease due to the shift of the 
Andreev states to higher energy. 
Phenomenological theory: Lifshitz invariants 

We now turn to a phenomenological description of the 
above microscopic results. This phenomenological theory 
shows that the appearance of a spin-triplet component is 
generic and not specific to the microscopic details. The 
key point is that the admixture of spin-singlet and spin- 
triplet order parameters is due to the existence of Lifshitz 
invariants in the Ginzburg Landau free energy (such in- 
variants where first discussed by Mineev and Samokhin 
25|). In particular, if the spin triplet order parameter 
has the form d{k,R) — J2i j ^iji-^)^^^] ^^'^ il i^s{R) 
describes the s-wave pairing, then symmetry allows the 
following LI (note that a similar LI has been found in the 



context of cold atoms |18| ) 



part becomes non-zero, these states are shifted away from 



J2 H,[A.,{i\/,yj,r + A* ,(*v,V'.)]. (12) 
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If, for example, the spin-singlet order parameter is given 
by V'(-R) = V'o cos{qRj), then this term implies Aij{R) = 
tptQiHiqsm{qRj) with ipto 0. This LI ensures that a 
spin-triplet component is always induced. This captures 
some of the main resuhs found in the microscopic the- 
ory: the triplet order parameter is largest where the spin- 
singlet order parameter vanishes; and the relative phase 
between the spin-singlet and spin-triplet order parame- 
ters is 7t/2. The LI can also be generalized to unconven- 
tional spin-singlet order parameters and the role of SOC. 
For example, if the spin-singlet pairing is d^2_y2, then 
the following Lifshitz invariant exists 



-A^,y{^Vy^pd)* 



(13) 



This implies /-wave spin-triplet pairing appears and once 
again, the magnitude of the /-wave component is largest 
where the d-wave component vanishes. This has been 



argued to be relevant in the organic (TMTSF)2X[17|,|26|. 
Furthermore, for example, in a tetragonal material with 
spin-singlet s-wave order ips , spin-orbit interactions allow 
the following LI 



(14) 



where tpi is defined through d{k,R) = {ipy{R)kx — 
'ipx{R)ky)z. In this case, the triplet component will have 
the same phase as the s-wave component (as opposed to 
the 7r/2 phase shift for the field induced LI). 

The existence of the LI plays another role not tied 
to the LO phase. In particular, it has been argued 
that a singlet to triplet phase transition may occur in 
(TMTSF)2X superconductors without the existence of 
a FFLO phase [s, [l3|- Such a transition is typically first 
order. The LI terms will transform this first order transi- 
tion into a pair of second order transitions between which 
lies an inhomogeneous singlet-triplet mixed phase. To 
understand this, consider adding the following simplified 
free energy to the LI in Eq. [T3] 



+/3p(|V.P + |Vyn' 



(15) 



For this free energy, without the LI, the singlet to triplet 
phase transition is first order. Near the normal to su- 
perconducting phase boundary, where it is sufficient to 
consider the quadratic free energy, the singlet to triplet 
transition will occur when as = at- Close to this point, 
when the LI is included, the quadratic free energy is al- 
ways minimized by introducing a inhomogeneous state 
where both ips oc e"^'' and -0^ cx e*''". This solution in- 
tervenes between the pure singlet and triplet states and 



the transition into this inhomogeneous phase is second 
order from both the pure singlet and pure triplet phases. 
This indicates that even if there is no FFLO phase in 
(TMTSF)2X, a closely related inhomogeneous singlet- 
triplet phase is likely to appear. 

In conclusion, we present microscopic arguments that 
show that the spatial line nodes of spin-singlet LO phases 
are removed by the appearance of a spin-triplet com- 
ponents. We show that this can be understood phe- 
nomenologically through the existence of Lifshitz invari- 
ants in the free energy which also ensure that the spin- 
triplet component always appears in a spin-singlet FFLO 
phase. This or related inhomogeneous singlet-triplet 
mixed states are likely to exist in the organic supercon- 
ductors (TMTSF)2X. This work is supported by NSF 
grant DMR-0906655. 
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